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Abstract. In this paper we introduce a notion of viscosity solutions for Eikonal equa- 
tions defined on topological networks. Existence of a solution for the Dirichlet prob- 
lem is obtained via representation formulas involving a distance function associated to 
the Hamiltonian. A comparison theorem based on Ishii's classical argument yields the 
uniqueness of the solution. 



1. Introduction 

Several phenomena in physics, chemistry and biology, described by interaction of differ- 
ent media, can be translated into mathematical problems involving differential equations 
which are not defined on connected manifolds as usual, but instead on so-called ramified 
spaces. The latter can be roughly visualized as a collection of different manifolds of the 
same dimension (branches) with certain parts of their boundaries identified (ramification 
space). The simplest examples of ramified spaces are topological networks, which basi- 
cally are graphs embedded in Euclidean space. The interaction among the collection of 
differential equations describing the behavior of physical quantities on the branches is de- 
scribed by certain transition conditions governing the interaction of the quantities across 
the ramification spaces. 

From a mathematical point of view, the concept of ramified spaces has originally been 
introduced by Lumer [19] and has later been refined and specified by various authors, 
e. g., J. von Below and S. Nicaise [20]. Since 1980, many results have been published 
treating different kinds of interaction problems involving linear and quasi-linear differen- 
tial equations (confer for instance Lagnese and Leugering [15], Lagnese, Leugering, and 
Schmidt [16], von Below and Nicaise [6]). 

As far as we know, fully nonlinear equations such as Hamilton- Jacobi equations have 
not yet been examined to a similar extent on ramified spaces. In the present paper we 
attack the problem by extending the theory of viscosity solutions for Hamilton- Jacobi 
equations to topological networks. The major task in this context is to establish the 
correct transition conditions the viscosity solutions are subjected to at transition vertices. 
As a matter of fact, these transition conditions make up the core of our theory, as they 
constitute the major difference from the classical theory of viscosity solutions. 

The main result of the present paper consists in the observation that the concept 
of viscosity solutions can indeed be appropriately extended to the class of first order 
Hamilton- Jacobi equations of Eikonal type on a topological network. Any generalization 
of existing concepts to new scenarios has to be justified by the preservation of essential 
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features. In the case of the theory of viscosity solutions, these features are uniqueness, 
existence, and stability We will show that our generalization of viscosity solutions to 
networks will be just as "weak" to yield existence, while being sufficiently "selective" in 
order to ensure uniqueness and stability with respect to uniform convergence. We will also 
demonstrate that our definition arises as a natural selection principle, which in particular 
selects the distance function as the unique viscosity solution of the Dirichlet problem for 
the standard Eikonal equation on networks. 

A different attempt to study Hamilton- Jacobi equations on networks has already been 
made in [1]. However, the aim of this paper deviates from the one addressed in the present 
paper: its main issue is to characterize the value function of controlled dynamics in R 2 
restricted to a network. Therefore, the choice of the Hamiltonian, which may be discon- 
tinuous with respect to the state variable, has to be restricted by assumptions ensuring 
both a suitable continuity property with respect to the state variable and the fact that 
the set of admissible controls be not empty at any point of the network. Additionally, 
the definition of viscosity solution characterizing the value function is different from our 
approach, as it involves directional derivatives of test functions in M 2 along the edges. 
In the present paper, Hamilton- Jacobi equations and differentiation along the edges are 
given in an intrinsic way making use of the maps embedding the network in M. N , hence the 
approach is intrinsically 1-dimensional. Moreover in our approach appropriate assump- 
tions at the transition vertices guarantee the continuity of the Hamiltonian with respect 
to the state variable. 

The existence of a viscosity solution is obtained by a representation formula involving 
a distance associated to the Hamiltonian (see [7], [11], [13] for corresponding results 
on connected domains), the solution turning out to be the maximal subsolution of the 
problem. Uniqueness, on the other hand, relies on a comparison principle inspired by 
Ishii's classical argument for Eikonal equations [14]. In this respect, the existence of a 
strict subsolution plays a key role. 

An important and classical problem in graph theory is the shortest path problem, i.e. 
the problem of computing in a weighted graph the distance of the vertices from a given 
target vertex ([4]). The weights represent the cost of running through the edges. A mo- 
tivation of our work is to generalize the previous problem to the case of a running cost 
which varies in a continuous way along the edges. In this case the aim is to compute the 
distance of any point of the graph from a given target set and this in practice corresponds 
to solve the Eikonal equation \Du\ = a(x) on the network with a zero-boundary condition 
on the target vertices. Moreover Hamilton- Jacobi equations of Eikonal type are impor- 
tant in several fields, for example geometric optics [5], homogenization [8, 18], singular 
perturbation [2], weak KAM theory [9, 10, 11], large-time behavior [12], and mean field 
games theory [17]. 

The paper is organized as follows: In section 2 we introduce the definitions of topological 
networks and viscosity solutions. In section 3 we collect some basic properties of viscosity 
solutions, in particular stability with respect to uniform convergence. Section 4 is devoted 
to the study of a distance function associated to the Hamilton- Jacobi equation, while 
section 5 presents the proof of a comparison principle. In section 6 the representation 
formula for the solution of the Dirichlet problem is given. 
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2. Assumptions and preliminary definitions 

We start with the definition of a topological network. 

Definition 2.1. Let V = {vi, i E 1} be a finite collection of pairwise different points in 
R N and let j E J} be a finite collection of differentiate, non self-intersecting curves 
in M. N given by 

TTj : [0,lj] ->■ R N , lj > 0, j G J. 
Set tj := 7Tj((0, lj)), Ej := 7r,-([0, lj]), and E :— {ej : j E J}. Furthermore assume that 

i) TTj(0),TTj(lj) G V for all j G J, 

ii) #{<s. n V) = 2 for all j G J, 

iii) Ej fl efc C V , and j^{ej fl e^) < 1 for all j, k E J, j ^ k. 

iv) For all v, w E V there is a path with end points v and w (i.e. a sequence of edges 
{ej}^ =1 such that #(Ej fl Ej + ±) = 1 and v G E\, w G e^). 

Then T := IJjej^j c ^ N ^ s called a (finite) topological network in R N . 

If V{ G Vdej we say that Cj is incident to Vi (ej incvj in short). For i G / we set Jncj := 
{ j G J : inc fj}. Observe that the parametrization of the arcs ej induces an orientation 
on the edges, which can be expressed by the signed incidence matrix A = {aij} it j € j with 

{1 if Vi G Ej and 7r,-(0) = 
— 1 if Vi G Ej and 7r,-(Zj) = Uj, 
otherwise. 

In the following we will study boundary value problems on f . Given a nonempty set 
Jb C /, we define dT := {vi, i G Is} to be the set of boundaries vertices, while for 
ir := I\Ib we call {vj, i G It} the set of transition vertices. We also set T := f \ <9r. 
We always assume i E Ib whenever #(lnci) = 1 for some i E I. We remark that in 
applications such as the shortest path problem it is interesting to impose a boundary 
condition also at the internal vertices, i.e. vertices with #(JnCj) > 1. 

Consider the subspace topology induced to f by M. N . It coincides with the topology 
induced by the path distance 

(2.2) %,x):=inf(/" ^(s^ds : t > 0, 7 G 1 for x, 7/ G f, where 



wo j 

i) 7 : [0, t] — > T is a piecewise differentiable path in the sense that there are t := < 
ti < • • • < t n+i := i such that for any m — 0, . . . , n, we have 7([t m , i m +i]) C e Jm 
for some j m G J, tt^ 1 o 7 G C 1 ^, f m+i ), and 

ii) j;E is the set of all such paths with 7(0) = y, *y(t) = x. 

For any function u : f — > H. and each j G J we denote by ■u J the restriction of u to e"j, i.e. 



it J := u 



o 7Tj : [0, lj] ->■ R. 
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We say that u is continuous in T and write u G C(f) if u is continuous with respect to 
the subspace topology of f . This means that v? G C([0, lj] for any j G J and 

M- 7 (7r~ 1 (i'i)) = u fc (7r j ^ 1 (i;j)) for any z G /, j,k G JnQ. 

In a similar way we define the space of upper semicontinuous functions USC(F) and the 
space of lower semicontinuous functions LSC(f ), respectively, and the space C(T). 
We define differentiation along an edge Cj by 

Q 

dju(x) := djU^itj 1 ^)) = —u-'(nj 1 (x)), for all x G 6j, 
and at a vertex vi by 

Q 

dju{vi) := djU 3 (ttT 1 (vi)) = ^^(^j 1 ^)) for 3 € Inc i- 

A Hamiltonian H : f x R — > R of eikonal type is a collection (H j )j e j with if J : [0, lj] x R — > 
R satisfying the following conditions: 

(2.3) FeC ([0,yxl), jGJ, 

(2.4) H 3 (x,p) is convex in p G R for any x G [0, /j], j G J, 

(2.5) H J (x,p) — >■ +oo as |p| — > oo for any x G [0, j G J, 

(2.6) Hi{Trj\vi),p) = H k (Ti k l (vi),p) for any p G R, i G J, j, G /nq, 

(2.7) H^n-^Vi)^) = H^Tr-^Vi), -p) for any p G R, i G /, j G JriCj. 

Remark 2.1. Assumptions (2.3)-(2.5) provide standard conditions in the theory of viscos- 
ity solutions (see [11], [13],). Assumptions (2.6) -(2.7) represent reasonable compatibility 
conditions of H at the vertices of T , i.e. continuity at the vertices and independence 
of the orientation of the incident arc, respectively (the network is not oriented). A typ- 
ical example of a Hamiltonian satisfying these assumptions is H^(x,p) := p 2 — a(x), 
j G J, where a(x) = a^{j(j l {x)) for x G e~j and a 3 G C°([0, lj]), a^(x) > for x G f, 
a 3 (ttJ 1 (vi)) = a k (n^ 1 (vi)) for any i G /, j,k G Inci. 

Definition 2.2. Let ip G C(T). 

i) Let x G 6j, j G J. We say that ip is differentiate at x, if <pP is differentiate at 
nj\x). 

ii) Let x = Vi, i G It, j, k G Inc^, j ^ k. We say that ip is (j, k)- differentiate at x, if 

(2.8) (i, J i)jp J [z j 1 (./•)) + a ik d k ip k (-K k 1 {x)) = 0, 

where (a^) as in (2.1). Moreover, we say that ip is differentiate at x if ip is 
(j, k)- differentiate at x for any j, k G 7nc„ j ^ k. 

Remark 2.2. Condition (2.8) demands that the derivatives in the direction of the incident 
edges j and k at the vertex Vi coincide, taking into account the orientation of the edges. 

On topological networks we now introduce the definition of viscosity solutions u of 
Hamilton- Jacobi equations of eikonal type of the form 

(2.9) H(x, Du) = 0, xeT. 
Definition 2.3. 

A function u G USC(T) is called a (viscosity) subsolution of (2.9) in T if the following 
holds: 
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i) For any x G ej, j G J, and for any ip G C(T) which is differentiable at x and for 
which u — <p attains a local maximum at x, we have 

W( 7 r- 1 (x),d m (7:- 1 (x))<0. 

ii) For any x = v^, i G It, and for any ip which is (j, k)- differentiable at x and for 
which u — <p attains a local maximum at x, we have 

H\nr\ x ),d m (nr\ x ))<0. 

A function u G LSC(T) is called a (viscosity) supersolution of (2.9) in T if the following 
holds: 

i) For any x G ej , j G J , and for any <p G C(T) which is differentiable at x and for 
which u — <p attains a local minimum at x, we have 

W{ixj\x),d m {*j\x))>0. 

ii) For any x = Vi, i G It, j G InCi, there exists k G InCi, k ^ j , (which we will call 
i-feasible for j at x) such that for any <p G C(T) which is (j, k)- differentiable at x 
and for which u — ip attains a local maximum at x, we have 

W(nr\ x ),d m (nr\x))>0. 

A continuous function u G C(T) is called a (viscosity) solution of (2.9) if it is both a 
viscosity subsolution and a viscosity supersolution. 

Remark 2.3. i) Leti G It andip G C(T) be (j, k)- differentiable atx. Then by (2.6)-(2.7) ; 

we have 

(2.10) W(n- 1 (x),dj l p(n- 1 (x))) = H k (n k \x), ±d k <pfa x (x))); 

hence in the subsolution condition, it is indifferent to require the condition for j or for k. 

ii) To simplify the notation, we set 

Hi(x,djip(x)) := W(7T- 1 (x) : d j <p(7r- 1 (x))). 

Moreover, we will call ip G C(T) an upper (lower) (j, k)-test function of u at x = if it 
is (j, k)- differentiable at x and if u — ip attains a local maximum (minimum) at x. 

iii) It is important to observe the asymmetry in definition 2.3 regarding the subsolution and 
the supersolution conditions at the transition vertices. It reflects the idea that distance 
functions have to be solutions of (2.9) and that there is always a shortest path from a 
transition vertex to the boundary. In fact it is worthwhile to observe that if supers olutions 
were defined similarly to subsolutions, the conditions in general would not be satisfied by 
m.t{d{y,x) : y G <9r} ; which is, as we will see in section 6, the solution of |Dm| 2 — 1 = 
with zero boundary conditions. 

iv) Taking (2.10) into account, it is easily seen that a viscosity solution u of (2.9) satisfies 
the equation in a pointwise sense at any point x G T where it is differentiable. 

As the conditions in definition 2.3 are of pointwise character, they can also be imposed 
on subsets of V. Hence, for Q C T we denote by S(Q) (S + (Q) or S~(Q), respectively) the 
space of solutions (supersolutions or subsolutions, respectively) of (2.9)in Q. For short, 
we set S := S(T), S + := S + (T), and 5" := «S~(r). 
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3. Some basic properties of viscosity solutions 

In this section we discuss some basic properties of the viscosity solutions introduced in 
the previous section. 

Proposition 3.1. Let u, v be subsolutions (supersolutions) of (2.9) in V. Then w := 
max{«, 11} (w := min{u,v}) is a subsolution (a supers olution) of (2.9) in T. 

Proof. We only consider the CclSG X — V i , % £ I) clS otherwise the argument is standard. 
Let u,d6 USC(r) be two subsolutions at x and observe that w = max{u,v} G USC(r). 
Let j, k G InQ and let ip be an upper (j, A;)-test function of w at x. If w(x) = u(x) 
(similarly in the other case) then ip is an upper (j, A;)-test function of u at x, implying 

H j {x,d jV {x)) < 0. 

Hence w is a subsolution. 

Let u, v G LSC(r) be two supersolutions at x, whence w = min{u, v} G LSC(T). Assume 
w(x) = u(x) (similarly in the other case). Hence for any j G Inci there exists k G InCi, 
k 7^ j, such that for any lower (j, A;)-test function <p of u at x we have 

H j (x,d jV (x)) > 0. 

Hence k is i-feasible for j also with respect to w. □ 

Proposition 3.2. Assume H n (x , p) — > H(x,p) uniformly for n — > oo (i.e. HiiixJ 1 ^)^) — > 
Hi {jxj 1 (x) , p) uniformly for (x,p) G e,j x M for any j G J). For any n G N let u n be a 
solution of 

(3.1) H n (x, Du) = 0, xeT, 

and assume uniformly in T for n — > oo. T/ien u is a solution of (2.9). 

Proof. We treat the case x = Vi, i G It, as otherwise the argument is standard (see [3]). 
We first prove that u is a subsolution. Choose any j, k G InCi, j ^ k, along with an upper 
(j, /c)-test function ip of u at x. Consider the auxiliary function (ps(y) '■= <p(y) + Sd(x,y) 2 
for S > 0. Observe that d m (d(x, •) 2 )( 7r m 1 ( ;r )) = for m = j and m = k, hence d(x, -) 2 
is (j, /c) differentiable at x. Then <ps is an upper (j, /c)-test function of w at x and there 
exists r > such that u — ip$ attains a strict local maximum w.r.t. B r (x) at x, where 
B r (x) := {y G T : d(x, y) < r}. Observe that a; is a strict maximum point for u — <ps also 
in £> := S r (x) fl (e^ U e^). Now choose a sequence u n — > for n — )• oo with 

(3.2) SUp |-u(x) — Mri(x)| < UJ n 

r 

and let y n be a maximum point for u n — ips in B. Up to a subsequence, y n ^ z E B. 
Moreover, 

m(x) - ^(x) - u n < u n (x) - tp s (x) < u n {y n ) - (p 5 (y n ) < u(y n ) - ip s (y n ) + u n . 

For n — )> oo, we get u(x) — <ps(x) < u{z) — (ps(z). As x is a strict maximum point, we 
conclude x — z. Invoking 

u(x) + (p s (y n ) - ips(x) -uj n < u n (y n ) < u(y n ) + u n 

we altogether get 

(3.3) lim y n = x, lim u n (y n ) = u(x) 

n— >oo n— >oo 
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We distinguish two cases: 

Case 1: y n ^ x. Then y n G e m with either m = j or m = k. Since u n — ips attains a 
maximum at y n and as (p$ is differentiable at y n with d m <ps(yn) — d m Lp(y n ) + 2Sa im d(x, y n ), 
we have 

(3.4) H™(y n , d m <p(y n ) + 25a tm d(x, y n )) < 0. 
with either m = j or m = k. 

Case 2: y n = x. Then d m ips(y n ) = d m <p(y n ) for m — j and m = k and therefore 

(3.5) Hl(y n ,d m (y n ))<0. 

By (3.3), (3.4), (3.5) and recalling (2.10) we get for n — > oo 

H j {x,d jV {x)) < 0. 



To show that u is a supersolution, we assume by contradiction that there exists j G IriCi 
such that for any k G Inci, k ^ j, there exists a lower (j, fc)-test function ipk of u at x for 
which 

(3.6) H*(x,dj<p k (x))<0. 

By adding a quadratic function of the form —akd(x, y) 2 to the function ip k we may assume 
that there exists r > such that u — ip k attains a strict minimum in B r (x) at x. Observe 
that x is a strict minimum point of u — ipk also in B k := B r {x) fl (ej Uej). Now for any 
nGN there exists A; n G iraQ, A; n ^ j, which is 2-feasible for j with respect to u n . Up to a 
subsequence, we may assume that there exists k G Inci such that k n = k for any n. 
Let y n be a minimum point of u n — ipk in B k and let u n be as in (3.2). Similarly to the 
subsolution case, we can prove that (3.3) holds. If y n ^ x, we obtain 

(3.7) H™(y n ,d m (yn))>0 
for either m = j or m = k. If y n — x, we get 

(3.8) Hi(y n ,d jiPk (y n ))>0. 
Hence by (3.3), (3.7), and (3.8), we get for n — > oo 

H^x,djip k ^))>0, 
which is a contradiction to (3.6). □ 

The proof of the following proposition is given in [3, Prop.II.4.1], for example. 

Proposition 3.3. Let K be a compact subset ofT and let u G S~(K). Then there exists 
a constant Ck depending only on K such that 

(3.9) \u{x)-u{y)\<C K d{x,y). 

The proof of the next two propositions is very similar to the one of Prop. 3. 2. 

Proposition 3.4. Let T C S~ (T C S + ) and set u(x) := sup{v(x)\v G T} (u{x) : = 
mi{v(x)\v ET}) for x G T. Suppose that u G C(Y). Then u G S~ (u G S + ). 
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Proof. To prove that u(x) = sup{v(x) \ v G T} is a subsolution, we only consider the case 
x = Vi, i G It- Consider j, k G IriCi, j ^ k, and an upper (j, fc)-test function ip of u at x. 
Set (ps(y) '■= f(y) + Sd(x, y) 2 for 5 > 0. Then ^ is an upper (j, /c)-test function of w at x 
and there exists r > such that u — ip$ has a strict local maximum point in B r (x) at x. 
Observe that x is a strict maximum point for u — <p& also in S := B r (x) fl (e^ U e^). Let 
u n G 5 be such that 

— < — 

n 

and let y n be a maximum point for u n — ^ in B. Up to a subsequence, y n — > 2; G 5. 
Moreover, 

u(a;) - - - < u n (x) - ^(x) < u n (y n ) - (ps(y n ) < u(y n ) - (f S (y„). 

n 

For n — > 00, we obtain — tps(x) < u(z) — <fs( z ), implying x = z, as x is a strict 
maximum point. Moreover, by 

u(a;) + <ps(y n ) ~ <Ps(x) < u n (y n ) < u(y n ) 

we get 

lim y n = x, lim u n (y n ) = u(x), 

n-^-oc n— s-oo 

and we conclude as in Proposition 3.2. 

Similarly to the proof of Proposition 3.2 one can also show that u(x) := inf{v (x)\ v G T} 
is a supersolution. □ 

Proposition 3.5. Let T C S and let u(x) := inf{t> (x) \ v G T} /or x G T. Assume that 
u(x) G R /or some xGT. T/ien u E S. 

Proof. By (3.9) all r G <S are uniformly Lipschitz continuous. As u(x) G R, we thus have 
u(y) G R for any y G T. Moreover, u is Lipschitz continuous on T. Next observe that by 
Proposition 3.4 u is a supersolution of (2.9). 

In order to prove that u is also a subsolution we once more invoke (and only sketch) 
the argument used in Proposition 3.2. Consider x = v iy i G It, j, k G InCi, j ^ k, 
and an upper (j, fc)-test function ip of u at x. Define the auxiliary function (ps(y) '■ = 
<p(y) + 5d(x,y) 2 for 5 > 0. Then ^ is an upper (j, fc)-test function of u at a; and there 
exists r > such that u — Lp$ has a strict local maximum point in B r (x) at x. Observe 
that a; is a strict maximum point for u — tpg also in B := B r (x) fl (ej U e^). Let u n e T 
be such that 

-u(x) — M„(x) > 

and let y n be a maximum point for u n — ipg in S. Up to a subsequence, y n — >■ z G 5. 
Moreover, 

1 

- ^(x) < u n (x) - ifs(x) < u n (y n ) - (p s (y n ) < u(y n ) - (ps(y n ) + -. 
Hence we obtain (3.3). Arguing as in Proposition 3.2 we conclude H j (x, dj(p(x)) < 0. □ 
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4. A DISTANCE FUNCTION FOR HAMILTON- JACOBI EQUATIONS 

In this section we assume 

(4.1) 5"(r)^0, 

i.e. there exists a subsolution of (2.9) in T. We introduce a distance function related to 
the Hamiltonian H on the network. For x, y G T define 

(4.2) S(y,x) = My\( 1 (s)^(s))ds: t > 7 G , 

where 5* as in (2.2) and 

L J (x, q) := sup{pg — W(x,p)} = sup{pq — H 3 (ttJ 1 (x) , p)} 

peK PGR 

for any j G J, x G e,-. Note that the distance defined by (4.2) coincides with the distance 
defined by (2.2) for H(x,p) = \p\ 2 — 1. The next proposition summarizes some properties 
of S. 

Proposition 4.1. S is a Lipschitz continuous distance onT x T. Moreover, 

i) for anyyeT we have S{y, •) G S"(r) n S(T \ {y}), 

ii) for any x,y G T we have 

(4.3) S{y,x) = max{u{x)\ueS-{T)s.t. u(y) = 0}. 

Proof. By (3.9) any subsolution u of (2.9) in T is Lipschitz continuous. Integrating along 
a path joining x and y we get 

(4.4) u(x) — u(y) < S(y, x) for any x, y G T. 

Thus by (4.1) S^y, x) is finite for any x, y G f . 

By the coercitivity of H assumed in (2.5) there exists constant R > 0, M such that 
L(x,q) < M (i.e L 3 (nj 1 (x) , q) < M) for any x G T, q G B(0,R), and therefore S(y,x) < 
Cnd(y,x) (see [12, Prop. 5.1] for details). Moreover, given x,y,z G T, the juxtaposition 
of two curves in and B zx gives a curve in B yx , whence 

S(y,x)<S(y,z) + S(z,x). 

S(y, •) is a subsolution: In order to prove that S(y, •) is a subsolution at Xo G T we 
distinguish two cases: 

Case 1: Xo G" {vi,i G It}- Assume Xo G ej for some j G J and let -0 be an upper 
test function of S(y, •) at x . It follows S^y, x ) — ipj(nJ 1 (xo)) > <S'(y, x) — ipj(7rj 1 (x)) for 
x G B r (xo) fl ej. Set to : = Tr/ 1 ^), ^ x 9 e ^> an d choose ft. sufficiently small in such a 
way that t — hq <E (0, /j). Define the curve 7^ : [0, h] — > T by 7h(s) := 7T,-(| to + (1 — f )th), 
where th := t — hq and set Xh q '■= Kj(to — hq). Hence 

Ft K \ a 1 u \ r V^o) -Tpjjth) ^ v S{y,x ) - S{y,x hq ) 
dip(x ) q = djipjfo) q = km — : — < km 



< lim S(x h x ) 1 t L{ {s) %{s))ds = 

h^0+ h h^0+ h J 



1 f h . s s 
km - / L J (-t + (1 - -)t h ,q)ds = L 3 (t ,q) = L(x ,q). 
h^o+ h J h n 
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Hence H(x, dip(x )) = sup q {dip(x ) q — L(x , q)} < 0. 

Case 2: x G {vi,i G It}- Assume x = Vi and let ip be an upper (j, /c)-test function 
of S(y, ■) at xo. Set t m := n^ m 1 (xo), m = j, k, and observe that for any q G R and for h 
sufficiently small we have t m —h(—ai m q) G (0, l m ) (or equivalently 7r m (t m — h(— cii m q)) G e m ) 
for m — j, k. Arguing as in case 1, we get 

dj^j^ixo^i-ciijq) < LP(t:J\xq), -o^g), 

dk^k{^ k (x ))(-a ik q) < L (ix k (xo),-a ik q). 

Moreover, since a ij d j ii(n~ 1 (xQ))+a ik d k ^ k ('K k ~ 1 (x )) = and L J (ttJ 1 (x ) , q) = L k (n^ 1 (x ) , q) , 
we get 

^(^(xo))^) < L^tt -1 ^),^), 
d k ^ k (ir k (x ))(a ik q) < L (n k (x ),a ik q). 

By (4.5) and (4.6) it follows d m ip m ('K^(xQ))q < L m '((^m ( x o) > Q) f° r m = 3\k, whence 
H(x, dip(x )) = sup {dtfj(x )q — L(x , q)} < 0. By (4.4) and since S(y, •) is a subsolution 
with S(y,y) = 0, we finally obtain (4.3). 

S(y, •) is a supersolution inT \ {y}: In order to prove that S(y, •) is a supersolution at 
x 7^ y we only consider the case x — Vi, i G It, being the other case standard. Assume 
that u(-) = S(y, •) is not a supersolution at x. By definition there exists an index j G Inci 
for which there does not exist any i-feasible index k G InCi, k ^ j. Hence for any 
k G K = IriCi \ {j} there exists a lower (j, fc)-test function ip k of u at x with 

(4.7) W{x,dj<p k {x)) <0. 
By Remark 2.3(i) we have 

(4.8) H*{x, djip k {x)) = H k (x, d k ip k (x)) < 0. 

It is not restrictive to assume that u(x) = <p k (x) for any k G K. By adding a term of the 
form — ad(z, x) 2 we may assume that u — ip k attains a strict minimum point at x. Hence 
by (2.3) and (4.8) there exists r > such that for all k G K and z G B r (x) \ {x} 

u(z) - tpk(z) > 

(4 ' 9) H j (z, d m {z)) < 0, H k (z, d k ip k {z)) < 0. 

Let £ > be such that 

(4.10) u(z) - ip k (z) > £ for all k G K and z G dB r (x). 

Define ^(z) := <^(z) + £ and v : {x} U \J k£lnCi e k -> R by 



max^^(z), if z G e^, 
(p k (z), if z G e fc , k G X. 



We claim that i> is a subsolution of (2.9) in B r (x). 

Case 1: Consider z G B r (x) fl e\. If / G then v(z) = <fii(z) and the claim follows by 
(4.9). If I = j, by (4.9) we have 

H k {z,d k <p k {z)) <0 V&gK 
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and the subsolution condition follows by Proposition 3.1. 

Case 2: Consider z — x. First assume l,m G K, I ^ m, and let if) be an upper 
(I, m)-test function of v at x. Set 

di := audiipix), d rn := a im d m ip(x), rji := a u di(fii(x), r} m := a im d m <p m (x). 

As if) is (I, m)-differentiable at x, we have di + d m = 0. If di < 0, we have d\ > rji by 
the definition of v and by the fact that v — if) attains a local maximum at z — x. Hence 
\di\ < \r)i\. Similarly, if d m < 0, we have d m > r) m , implying \d m \ < \r) m \. We therefore 
conclude that 

1^0*01 = \9mip(x)\ < max{|^^(x)|, |<9 m <^ m (a;)|}. 

By the assumptions on H (see (2.3)-(2.7)) the function h : R — >■ R, p i— > H s (x,p) for 
s G Inci is independent of s, symmetric at p — 0, and strictly increasing in Hence by 
(4.8) 

H l (x,ditp(x)) = h(diif)(x)) < max{h(di<pi(x)), h(d m <p m (x))} = 
max{H l (x,di<Pi(x)), H m (x,d m fi m (x))} < 0. 
Assume now that if) is an upper (J, /)-test function of v at x and set 

dj := aijdiif)(x), d t := a a diif)(x), ej := max aijdjip k (x), rji := a a di(pi{x). 

As above, we have \dj\ < \ej\ if dj < and \d m \ < \rj m \ if d m < 0. Hence we get 
\djif)(x)\ = \diif)(x)\ < max{max\d k (p k (x)\, \di<pi(x)\}, 

and therefore by (4.8) 

H^{x,djif){x)) = h(djif)(x)) < max{maxh(dj<f k (x)), h(di<pi(x))} = 

k<=K 

max{maxH k (x,djif>i(x)), H l (x,di(fi(x))} < 0. 

Hence v is a viscosity subsolution in B r (y). Define the function v : T — > R by 

q ,( y \ — / ™x{v(z),u(z)}, if z G Bt(x), 
{ > ■ \ u(z), if z G T\B t (x). 

By (4.10), v is continuous, v = u outside B r (x), and v is a subsolution of (2.9) in T. Since 
v(x) = v(x) > u(x), we get a contradiction to (4.3). □ 

5. A comparison Theorem 
This section is devoted to the proof of a comparison theorem for problem (2.9). 

Theorem 5.1. Assume that there exists a closed subset K C V and a function f G C(T) 
with f(x) < for all x G T \ K. Moreover, let u be a subsolution of 

(5.1) H(x,Du) = f(x), x eT\K, 

and let v be a supersolution of (2.9) inT\K . If u < v on dF U K , then u < v inT. 
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Proof. Assume by contradiction that there exists z G T \ K such that 

(5.2) u(z) — viz) = maxiw — v} = 5 > 0. 

r 

For e > define $ e : f x f ->■ R by 

$ £ (x, y) := u(x) - v(y) - e~ l d(x, y) 2 . 

As $ e is upper semicontinuous there exists a maximum point (p e , g e ) for $ e in f 2 . By 
$ e (z,z) < $ e (p e ,g e ) we get 

(5.3) t~ l d{p t , q e f < u{p e ) - v(q e ) - 5, 
whence 

(5.4) \imd(p e , q e ) = 0. 

By the compactness of T, there exists p G T such that — >■ p. By (5.3) and the 

Lipschitz continuity of u (see (3.9)) we get 

£~ L d(p t , qef < u{p e ) - u(q e ) + u{q e ) - v(q e ) - 5 < Ld(p e , q e ) 

and therefore 

(5.5) lim e~ 1 d{p e ,q e ) = 0. 

Moreover, by (5.3)-(5.4) we have p G T \ K as well as p t , q e G T \ K for a sufficiently small 
choice of e. Next observe that it is possible to assume that there exists a unique path 
7 = 7 e of length d(p e , q e ) in T connecting p e and q e which runs through at most one vertex 
Vi, i G I. We distinguish several cases (for simplicity we set p := p e , q := q e ). 

Case 1: There are indices i G / and j, k G IriCi such that p G e,-, g G e k and such 
that 7 runs through We observe that the functions ip p (x) := e~ 1 d(p, x) 2 and : = 

e^d(x, q) 2 are differentiable at q and p, respectively. In fact, if 

(5-6) P = nJ 1 (p), q = K k ~ l (q), 

we have 

d M(P) = ^d(p, q)a ij7 d k ipp(q) = e~ 1 d(p, q)a ik . 

Observe that u — ip q has a maximum point at p and v + ip p has a minimum point at q, 
whence 

Hi(p,d jVq (p)) = W(p,e- l d(p,q) ai] ) < f(p) 

H k (q,d k ip p (q)) = H\q,-e- l d{p,q)a ik ) > 0. 

We denote by u m , m = j,k, the modulus of continuity of H m with respect to (x,p) G 
e m x R. By (2.6) and (2.7) there is some rj > such that for sufficiently small e > we 
have 

V < ~f(p) < H k {q, -e^dip, q)a lk ) - H j (p, e~ l d{p, q)a ZJ ) < 
H k (vi, -e~ 1 d{p,q)a ik ) - H 3 (v i ,e~ 1 d(p,q)a ij ) + u k (d(vi,q)) +u j (d(v i ,p)) = 
H : '(v i ,^ l d(p,q)a ik ) - H 3 (v u e _1 d(p, q)a i:j ) + u k {d(v h q)) + u^v^p)) 
< u j (e' 1 d(p,q)(\a ij \ + \a ik \)) +u k {d{v u q)) +Uj{d{v h p)). 
By (5.5) we get a contradiction for e — > 0. 
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Case 2: There are indices i E I and j G Inci such that p G ej and q = V{. As q G T, we 
have i G It- Setting ip p and ip q as above and using a notation similar to (5.6) we have 

(5.7) djip? p (q) = -e -1 d(p, g)a^, <9 fe <^(g) = e^d(p, q)a ik for all G inq, fc 7^ j. 
Hence 

(5.8) dijdjippiq) + a ik d k (Pp{q) = (-a?- + a? fc )e -1 d(p, 5) = for all fc G /nq, fc 7^ j. 
Thus ip p is (j, &)-differentiable at q for all fc G iraq, k ^ j. Moreover, 

9 j^(p) = e -1 d(p,?)aij- 
Since u — ip q has a maximum point at p, if follows 

(5.9) /P(p, 9^,(p)) = if'(p, e- 1 d(p, gK) < /(p). 

Moreover, since t> +</? p has a minimum point at q — q, there is an i-feasible index k G InCi, 
k 7^ j, for j. By (5.8), y p is (j, /c )-differentiable, whence we obtain 

(5.10) H*(q, d jVp {q)) = W(q, -e" 1 ^, g)a lfe ) > 0. 

Subtracting (5.9) from (5.10) we derive a contradiction as in case 1. 

Case 3: There are indices i G / and j G InCj such that p = Vi and q £ ej. We proceed as 
in case 2, observing that the definition of subsolutions is less restrictive than the definition 
of supersolutions and therefore no extra argument is required. 

Case 4-' There are indices j G J such that p, q G and p 7^ g. Setting </? p and as 
above and using a notation similar to (5.6) we have 

(5.11) d^ q {p) = -d s <p> p (q). 
Hence we have 

Hl(p,d s <p q (p))<f(p), 

Hi(q,-d m (q)) = H^(q,d jVq (p)) > 

and we conclude as in the previous cases. 

Case 5: We finally assume that p — q- Assume p = q = Vi for i G It (the case p, q G ej 
for j G J is similar). Then 

a,^(vr- 1 K)) = a,^(7T7 1 (t; l )) = 

for all j G JnCj. In particular for each choice of j, k G InCi, both tp q and — <^ p are 
(j, /c)-differentiable and we get a contradiction as in the previous cases. □ 

6. Representation formula for viscosity solutions 
In this section we give a representation formula for the solution of the Dirichlet problem 

(6.1) H(x,Du) = 0, xeT, 

(6.2) u = g, x G dr. 

In addition to (2.3)-(2.7) and (4.1), in this section we assume that 

there exist a closed (possibly empty) subset K C T, a differentiable 

(6.3) function ip, and h G C(r) with /i < in T \ K such that 

H(x,Dip) < h{x), xeT\K, 
i.e. tp is differentiable the sense of Definition 2.2 and a strict subsolution in T \ K. 
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Proposition 6.1. Let g : T — > R be a continuous function satisfying 

(6.4) g(x) — g(y) < S(y, x) for any x, y G K U <9I\ 

where S is the distance defined in (4.2). Then the unique viscosity solution of (6.1) -(6.2) 
is given by 

u(x) := xmn{g(y) + S(y, x) : y G K U dT}. 

Proof By Proposition 3.5 u is a solution of (2.9). Observe that we have u(x) ^ g(x) for 
x G K U dT if and only if there is some z G K U dT such that g(x) > S{z,x) + g(z). 
However, this is ruled out by assumption (6.4). Hence u is a solution of (6.1)-(6.2). 
Assume that there exists another solution v of (6.1)-(6.2). For 9 G (0, 1) define Ug := 
9u + (1 — 9)ip, where ip as in (6.3). By adding a constant it is not restrictive to assume 
that ifj is sufficiently small in such a way that 

(6.5) ug(x) < u(x), x G T. 

First, let x G ej fl (r \ AT) for some j G J and let (/? be an upper test function of u at x. 
Setting ifo := 9ip + (1 — 6 1 )-?/; we obtain by means of convexity 

(6.6) H j (x, djipo) < 9H j (x, dj<p) + (1 - 9)H j (x, drf) < (1 - 6>)/i j (:r). 

Secondly, assume that x = Vi for some z G It- Fix any two indices j, k G InQ, j ^ k, 
and let be an upper (j, k)-test function of u at x. Setting tpg := 9tp + (1 — 9)ip and 
observing that by definition 2.2 tpg is an upper (j, A;)-test function of ue at x, we again 
obtain (6.6). Hence uq is a viscosity subsolution of 

H(x,dju) < (l-9)h j {x). 

Applying theorem 5.1 with / = (1 — 9)h and (6.5), it follows ug < v for all 9 G (0, 1). 
Letting 9 tend to 1 yields u < v. Exchanging the role of u and v we conclude that u = v 
in f . □ 

Remark 6.1. For the problem \Du\ 2 — a(x) = (see Remark 2.1) the existence of a strict 
subsolution follows by setting K := {x G T : a(x) = 0}, ip :— C for some suitable C G K, 
and h(x) := —a(x). 

If g does not satisfy assumption (6.4) we can still characterize S as the maximal solution 
of the problem. 

Proposition 6.2. Let g : T — >■ R be a continuous function. Then 

u(x) := mm{g(y) + S(x, y) : y <E K U dT} 

is the maximal solution of (6.1) among the solutions v of (6.1) which satisfy v < g on 
KUdT. 

Proof. By Proposition 3.5, u is a solution of (2.9). If v is a solution of (6.1), then by (4.4) 

v(x) < v(y) + S(y, x) < g(y) + S(y, x) for any y G K U dT, 
and therefore the statement follows by Theorem 5.1. □ 

Remark 6.2. As explained in the introduction, a motivation of our work comes from the 
shortest path problem on a network. The case of a weighted graph studied in graph theory 
fits in our framework. In fact it is sufficient to choose the function otj in Remark 2.1 in 
such a way that its integral along the edge ej is equal to the given weight. We will study 
this problem in more details in a forthcoming paper. 
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